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Abstract 

Wave propagation in curved tubular domains is considered. A 
general version of Webster's equation is derived from the scattering 
passive wave equation. More precisely, it is shown that planar aver- 
ages of a sufficiently smooth solution of the wave equation satisfy the 
corresponding Webster's equation when the latter includes additional 
control signals determined by the solution. 
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1 Introduction 

We study wave propagation in narrow and long, tubular 3D domains ft of 
finite length, with varying cross-section area. In this case, there is a classical 
approximation in one spatial variable, Webster's equation, to the full 3D 
wave equation. Variable cross sections of Q are represented there by the 
area function A{-) that appears as a weight in the elliptic part of Webster's 
equation. The purpose of this and the companion papers |3l |T5] is to justify 
this approximation in great generality and in a mathematically rigorous way, 
and to give an error estimate by an a posteriori method. 

Early work concerning Webster's equation can be found in [3, [251 EEl EE] . 
Webster's original work [28] was published in 1919, but the model itself has 
a longer history spanning over 200 years and starting from the works of 
D. Bernoulli, Euler, and Lagrange. A selection of more modern approaches 
is provided by [121 [131 [13 [211 [22], and [22], which is particularly interest- 
ing in the context of this work. Shared features of the modern treatments 
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are the time-harmonic Ansatz and asymptotic expansions for the velocity 
potential function satisfying the wave equation (either with or without 
underlying convection). The time-harmonic Ansatz leads to the correspond- 
ing time-invariant form of the model, analogous to the Helmholtz equation. 
The resonance structure (i.e., the formants) of the model is obtained from 
the resulting eigenvalue problem; see, e.g., asymptotic spectra of Neumann- 
Laplacian on shrinking tubular domains in [TTl [2l]. Our complementary 
approach is to look for equations that are satisfied by planar averages of 
which is likely how Webster's equation was first discovered. This results in 



a precise equation for the averages 0. If three load terms (3.11)-(3.13) are 
neglected from this equation, Webster's equation is obtained in a generalised 
form. 

Our interest in Webster's equation stems from the fact that it provides 
a simple and tractable models of the acoustics of the human vocal tract 
as it appears during vowel utterance. This equation (as well as the wave 
equation) can be used as a part of a computational physics model; see, e.g., 
[21 El El [HI [To] and the references therein. The other parts are a mechan- 
ical model for the glottis and an exterior space model outside the mouth; 
see [H [2] . Following these ideas and motivations, a boundary control action 
is used to represent the acoustic field produced by the glottis source. Fur- 
ther applications of Webster's equation include modelling of water waves in 
tapered channels, acoustic design of exhaust pipes and jet engines for control- 
ling noise, vibration, and performance as well as construction of instruments 
such as loudspeakers and horns. 

We represent a tubular domain i7 C as follows. The center line of the 
tube can be thought of as a smooth curve 7 (of unit length) parametrised by 
arc length s. We assume that the cross section perpendicular to the tangent 
of 7 at the point 7(5) is a circular disk r(s) with the radius R{s). The 
boundary of Q then consists of the ends of the tube, r(0) and r(l), and the 
wall F of the tube. 

Lossless acoustic wave propagation in fl can be modelled by the wave 
equation 

= c^A(f) in Qx 1+, 
g = ^andpf=p on F(0) x 1+, (1.1) 
= on F(l) and §^ = on F 

where u denotes the unit exterior normal vector, is the (perturbation) veloc- 
ity potential, and the control/observation functions p = p{r, t) and v = f (r, t) 
are the pressure and (perturbation) velocity at the controlled end F(0), re- 
spectively. The sound speed in and the density of the medium are denoted 



2 



by c and p, respectively. The Dirichlet condition on r(l) represents an open 
end, and the Neumann condition on F represents a completely reflecting sur- 
face. Other boundary conditions (such as the boundary dissipation condition 



(3.2)) can be used as well. 

The lossless Webster's equation 



V{t) 



A{s) ds {^('^) ds) 

S(o,t), P{t)- 



for (s,t) e [0, 1] X M+, 
(0,t) and V^(l,t) = for t G 



d± 
^ dt 



(1.2) 



is a computationally and mathematically simpler wave propagation model 
than (1.1). Here A{s) is the area of the cross section r(s), and the solu- 



tion ip is Webster's velocity potential. Webster's equation (1.2) serves as a 
good approximation of (1.1) for wave lengths long compared to R{s). More 
precisely, the solution ifj approximates the averages 
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A{s] 



for s e (0, 1) and t > 



1.3) 



r(s) 



of the velocity potential given by (1.1). 



Let us discuss the control/observation boundary conditions on r(0) in 
(1.2). When either one of the models (1.1) and (1.2) is used, the wave 



propagation is treated in impedance form as in |H Section 5] . Alternatively, 
the control and observation boundary condition may be defined so that the 
physical dimension (which is power per unit area) of both input and output 
coincide. We take the latter approach in this paper as well as in [3l [15]. It is 
called the scattering form, and it is obtained by replacing the second line in 



(1.1) by the equations 



and 



;i-4) 



for t G M"*" and r G r(0) where u{r,t) is the input, and y{r,t) is the output. 
The interplay between impedance and scattering forms is discussed in a gen- 
eral setting in, e.g., [121 [13 [I8j. The corresponding inputs and outputs for 



the Webster's equation are uniquely obtained from (1.4) by averaging in the 



same way as in (1.3). Indeed, we get (under a minor technical assumption 
explained at the end of Section lil) 



-cf,0.,.|,0., 



dip 
~dt 



(0,t) 



and 



1.5) 
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for t e 



where u{t) 



AM/r(o)«('^'^)^^ and y{t) := ^ 
are the averages of the input u and output y, respectively. 

The wave equation with the boundary conditions = on r(l), equa- 



tions (1.4) and (3.2) define an internally well-posed, dissipative dynamical 
system as shown in [3l Theorem 4.1, items (i) and (iii), and Corollary 4.2]. 
The analogous statement holds for Webster's equation with the boundary 
conditions = and (1.5); see [31 Theorem 5.1]. Such systems have 

unique solutions and ip for twice differentiable input signals and compatible 
initial conditions at t = 0, as a consequence of [3| Proposition 2.2]. Further, 
by [31 Theorem 4.1(ii) and Corollary 4.2], has the regularity properties 



and A0 G C(M+;L2(fi)) 



(1.6) 



which will be used throughout this article. Let us note that the method of 
asymptotic expansions also requires a nontrivial implicit regularity assump- 
tion so that the expansion exists; see [221 P- 1985]. 

Referring to Fig. [T| the main results of this and the companion papers 
[H [15] can be summarized as follows: 



i) By Theorem 3.2 and the solvability of the scattering wave equation 



model on Q, all the signals in Fig. [Tj are well-defined functions, and the 
average satisfies nonhomogeneous Webster's equation with the load 



term F + G + H as given by (3.11) - (3.13) 



The solution of the homogenous Webster's equation ip approximates 
the averages in (1.3), and the tracking error can be controlled by F, 
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G, and H with the aid of a slight refinement of the argument in [21 
Section 5] to include direct control to state, as shown in [15] . 

The load terms F, G, and H (and hence, the tracking error) can be a 
posteriori estimated by using the technical results from Section |4| as 
shown in [T5l. 



All these results are shown for a more general model than the scattering 
version of (1.1) and (1.2). The curvature of Q is taken into account 



m 



Webster's equation by introducing a correction factor to the speed of sound; 



see (4.12) and also [22l Section 6]. Further, boundary dissipation is allowed 



through the lateral wall F by the boundary condition (3.2). This leads to 



a dissipation term in Webster's equation (3.4). The full results require that 
the torsion of the center line 7 of vanishes, and that ^'(0) = k(0) 
where k is the curvature of 7. 
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Figure 1: Feedforward coupling describing the main results of this paper. The 
equations in the blocks are as they appear in the lossless case and without 
curvature. 



The paper is organized as follows: In Section [2} we discuss the details 
of the coordinate system for the problem described above. The main result, 
Theorem 3.2 is given without proof in Section |3} In Section [4], we present 
the formal derivation of Webster's equation with curvature and dissipation. 
Also the scattering boundary conditions are derived from those of the wave 



equation. Section ^ contains the rigorous proof of Theorem 3.2 



2 Global coordinates in a curved tube 

In this section, we give the necessary facts about the coordinate system we 
use to represent tubular domains as shown in Fig. |2| see also [22], Section 2.3]. 
Detailed arguments are given in appendices at the end of the paper. 

Let 7 : [0, 1] —7- be a smooth curve parameterized by arc length s. 
Define for s G [0,1] the unit vectors 

t'(s) 

t(s) := 7'(s), n(s) := and his) := t(s) x n(s). 

k{s) 

Here the curvature of 7 is defined by k{s) := ||7"(s)]| for s G [0, 1], and we 
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Figure 2: The Frenet frame of the planar centrehne for a tubular domain fi, 
represented by some of its intersection surfaces r(s) for s G [0, 1]. The wall 
F C dQ is not shown. 

assume that miusgjo^i] k{s) > 0. It then follows that 

t' = nn, n = —nt + rb, and b' = — rn (2-1) 

where the number r = t{s) for s G [0, 1] is called the torsion of 7. All 
these facts concerning the Frenet coordinates can be found in any standard 
reference on elementary differential geometry, e.g., |20]. 

By i, j, and k denote the unit vectors in the direction of cartesian coor- 
dinate axes. For any point r = xi + + 2;k G in a neighbourhood of 7 
we have 

xi + yj + zk = 7(s) + r cos 9 n(s) + r sin 6* b(s), (2.2) 

and hence we can use the triple (s, r, 9) G [0, 1] x [0, 00] x [0, 27r) as coordinates 
for r. We shall call these coordinates {s, r, 9) the tube coordinates. 

Using these coordinates, we can define the (open) interior of any tube as 

n ■= {(s, r, ^) : s G (0, 1), r G [0, R{s)), 9 G [0, 2n)} C (2.3) 

where the differentiable function R : [0, 1] (0, 00) gives the radius of the 
cross section at s G [0, 1]. We define 

fr ■.= {{s,R{s),9):se{0A),9e[0,27T)} and 

\t{s) ■.= {is,r,9):re[0,Ro),9E[0,27r)} for sG[0,1]. 

As before, we denote by A{s) the area of T{s) for s G [0, 1]. We call F the 
wall, and the circular plates F(0), F(l) are the ends of the tube fl. The 
boundary of Q satisfies dfl = F U F(0) U F(l). We make it as a standing 
assumption that 

r]{s) := R{s)k{s) < 1 for all s G [0, 1] (2.5) 
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in which case the tube does not fold onto itself, and the coordinate mapping 



associated to (2.2) is a diffeomorphism. The number ri(s) is called curvature 
ratio at s G [0, 1], and typically both //(s) ^ 1 and ri'{s) ^ 1. The curvature 
factor is defined by 

- ^ (2.6) 



1 — rK(s) cos 9 



Clearly S = 1 + rKcosO + 0{ri'^) as r/ — t- since |S — 1 — rKcos6\ < 7/^(1 — 
For the derivation of Webster's model, we need the following facts about 



the coordinate system given by (2.2). We present the elementary computa- 



tions required to obtain these facts in the appendix. 

(i) The volume differential in tube coordinates is 

dV = dx dy dz = — ds dr d6. (2-7) 

(ii) The surface element on the tube wall V is 

dS = Wdeds where W{s) := R{s)^R'{sY + ((q{s) - If (2.8) 

(iii) The gradient V can be written as 

V = t{s)Di + n{s)D2 + h{s)D^ (2.9) 
where the differential operators Di, i = 1,2,3, are given by 

^ ^ d sinO d ^ ^ d cos 6 d 

D2 ■= cosO— — , D2, := smfe'TT — h 



dr r do dr r 89 

(iv) Finally, the exterior normal derivative on F is given by 

u.V=g (-R'E (l-rl]+il- .R)l] . (2.11) 



W \ \ds 06 J dr 
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3 Approximation of the wave equation 

In this section we present how the solutions of the wave equation 



c A(f) on the domain 



X Q 



(3.1) 



are related to a generalized Webster's equation in the domain that has been 



described in (|2.3|). We assume that the tube wall F, defined by (|2.4|), allows 

(3.2) 



energy dissipation through it, and the energy loss is modelled by 
d(j) d(t) 



^ dt dv 



on 



X r where a > 0. 



Note that the Neumann boundary condition ^ = 0, i.e., a = in (3.2), 
describes a perfectly reflecting boundary leading to a lossless model. For the 
purpose of reprensenting the curvature of the tube, we define the sound speed 
correction factor by 



'1/2 



where 77(5) is given by (2.5). (3.3) 



As shown in Section |4| Webster's equation is then given by 

2TTaW{s) dip 1 d^ip 



A{s) ds \ ds 



A{s) dt c2S(s)2 dt^ 



f 



(3.4) 



where / = 0. Note that the boundary dissipation ( |3.2 ) gives rise to the 
second term in (3.4). 



On the end surface r(l) of tube Q we use the Dirichlet boundary condition 

= on M+ X r(l) (3.5) 
that represents an acoustic open end. The same argument that produces 



from (3.5). 



(3.4) from (3.1) - (3.2) also produces the boundary condition ■ip{l,t) = 



Recall that the wave equation (3.1) is controlled and observed from the 



boundary component r(0) representing one end of the tube f2. In this work. 



we use the scattering boundary conditions given already in (1.4): 



^^''^ = 2^y(r,t) 



u{r,t) and 



(3.6) 
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Treating these boundary conditions is somewhat more comphcated that (3.5 ) 



and we make two additional assumptions on to ensure that the surface r(0) 
is "geometrically reflection-free" : 

A'{0) = and k{0) = 0. (3.7) 

As is shown in Section|4| then (3.6) imply the scattering boundary conditions 

d'ip d'ip 

" ' ' ;u,rj = ^^^/^^u[l) ana 

(3. 



-cg(0.)-f(0.) 



pA{0) 



m 



for Webster's equation. Other kinds of boundary conditions, even nonlinear 
ones, on r(0) are possible, and their treatment follows the lines of Section |4| 
The details are left to an interesting reader. 

The results of this section require the unique solvability of the wave equa- 
tion part in Fig. [TJ 

Proposition 3.1. Let the input signal satisfy u G C^(]R+; L^(r(0))) and let 
00 G H^{fl) satisfy 



dt 



dv ir(o)ur 



GL2(r(0)ur). (3.9) 



Then there is a unique solution of the wave equation (3.1), satisfying reg 



ularity conditions (1.6), the first (input) boundary condition in (3.6), the 



boundary conditions (3.2), (3.5), and the compatible initial conditions 



0(-,O) = 00, 



dt 



.0) 



dt 



(3.10) 



The output, given by the second equation in ( 3.6[ ), satisfies y E C(]R+; L^(r(0))) 



This is a reformulation of [3l Theorem 4.1(ii) and Corollary 4.2]. When 
consistent signals and initial conditions are used in Fig. [T] for the Webster 
part, we prove unique solvability there as well but this is not needed in this 
paper. 

Recall the feedforward coupling in Fig. [T] where both the wave equation 
and Webster's equation systems are boundary controlled by a common (apart 
from averaging over r(0)) external signal. To produce an exact cancellation of 
outputs in Fig. [1} we need to directly control the state of Webster's equation 
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by the forcing functions F, G, and H 



F{s,t) :-- 
G{s,t) :-- 
H{s,t) :-- 



1 d 



A'(s) 



A{s) ds 
27raW{s) d 
A{s) dt 



1 

2^ 



27r 



27r 



27r 



1 



V{s) 



-V - ■V(PdA- 



+ r[s] 



d_ 
ds 



V{s) 



d(j) 



dA - 



A{s 
aW{s)7]{s 
A{s) 



(P{s,R{s),e)d9jj ; (3.11) 
s,R{s),e)de^ ; and (3.12) 
EA(pdA (3.13) 



T(s) 



2n 



a' 



s, R{s), 9) COS 9d9 



where the error function is defined by 



E{s,r, 



1 



1 



2rK{s) cos 9 + K{sy{r' cos^ 9 - R{sY/A). (3.14) 



S2 S(s)2 

Now we are ready to give the main resuh of this paper. 

Theorem 3.2. Assume that : IR+ x — )■ M «s a solution of the wave equa- 
tion (3.1) satisfying the regularity properties (1.6) and the boundary condi- 
tions (3.2) and (3.5). 

(i) The averaged solution (j) : IR+ x (0, 1) — )■ M, defined by 

1 



A(.) 



r{s) 



satisfies the regularity conditions 



d(j) 



eC\R+;L\0,l)) and ^ e C\R+; L\0,1)), 

OS 



(3.15) 



(ii) The forcing terms F = F{s,t), G = G{s,t), and H = H{s,t) (as de- 
fined m ( [3ll| - ( |3l3| ) ; satisfy F e C{R+; L\0, 1)) and G G G\R+; L^{0, 
If m addition, r(s) = for all s m ([3l3|), then H e Ci(M+; L'^{0, 1)). 



(iii) //t(s) = m (3.13), i/ien t/ie averaged solution (p is a weak solution of 



1 d 
~A{s)¥s 



A{s) 



d(f)\ 27caW{s) d(f) 



ds 



A{s) dt c2E(s)2 dt^ 



(3.16) 



F + G + H. 
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(iv) // (j) satisfies the control/ observation boundary conditions (3.6) and 
R'{0) = k{0) = r(0) = 0, then satisfies the corresponding bound- 
ary conditions (3.8). 



This theorem is proved in Section |5j Claim (iii) means plainly that 



d<i) dC 
ds ds 



c2S(s) df' 



C A{s)dsdt 



{F{s, t) + G{s, t) + H{s, t))C{s, t)A{s)dsdt 



27raW{s)^Cdsdt 



(3.17) 



for all test functions C e C^HO, 1) x (0,T)) and all T > 0. Claims 1^ and 
(iv) state that the feed-forward connection in Fig. [l]is well-defined, and the 
output on the right vanishes. 

All the solutions and signals are functions by claims (pi) and (pTl). If equa- 



tions (3.4) and (3.8 ) together with t) = (corresponding ( |3.5[ )) are inter- 
preted as a scattering passive linear dynamical system (i.e., a system node) 
Sw on the state space = H^^y{0, 1) x L'^{0, 1) as in [3l Theorem 5.1], it 
is, in particular, a well-posed linear system in the sense of [HI Definition 

2.7], [27j. Thus, the function x : = 

Xyj as defined in fITi Definition 3.8. 

associated to S,„- 



P at 



is the unique strong solution in 
and Theorem 4.3.1] of the dynamics 



m 



s,. 



xit) 
u{t) 

fit) 



for almost all t > 



where f := F + G + H] see, e.g., [151 Section 2]. 



4 From wave equation to Webster's equation 



In this section, we give formal computations leading from the wave equation 
to Webster's equation. 



Derivation of Webster's equation 



To derive Webster's equation (3.4), we first obtain a weak formulation of the 



wave equation (1.1) - (3.2) in terms of tube coordinates. Open portions of 



the tube Vt are defined by 

fi(so,Si) = {(s,r,^) : s e (so,Si),r G [0,i?(s)),^ G [0,27r)} 
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for < So < si < 1. We denote r(so, si) := F fl n(so,si). On the cross 
sectional surface r(so) we have v = — t(so), and similarly we have z/ = t(si) 
on r(si). Thus by ( |2.9 ) - (2.6) we get for the exterior normal derivatives 



and 



ld_ 



d d 



d d 



for {so,r,9) eT{so) 



for {s,,r,9)eT{s,). 



(4.1) 



(4.2) 



Clearly dA = rdrdO on both r(so) and r(si). 



We use the Divergence Theorem and (4.1) - (4.2), and obtain 

■~-^^(t)dV 



Q(so,si) 

r2(so,si) 
r2(so,si) 
r(so) 



V(t)dV + 
y \^] - VcpdV -a 
dA + 



r(s„,si) 



H dp 



dA + 



r(so)ur{si) 



(90 d(p 



r(so,si 
r{si) 



1^ 



dA 



dA 



V 



n{so,si) 



+ 



r(si) 



(9s 



■VcpdV-a 
(90 



r(^o) 



9s 



r(so,si) ^ 
dA ] - f r(si) 



1 d(j) 



dA 



r(si) 



de 



dA - r(so) 



r(^o) 



09 



dA 



where Jj 



r(so,si) H du 



-a J^^^^^s^^dA is implied by (3.2). On the other 



hand, using the wave equation (|1.1|) and the volume element (2.7), we get 

2iT pR{s) q2a. 



I 

Jn 



A(PdV 



O{so,si) 



1 



dV 



n(so,si) 

Sl 



so 



1 



Sl 




Jo 



r dr dd 



V{s) 



dA ds. 



Combining these two expressions for the same thing, we get 

1 



L(so,Si) 



'n{so,si) 
+ r(si) 



V \^:^ ] -Vc^dV 
(90 



r(^i) 



dO 



dA — r(so) 



r(so) 



90 

de 



dA 



(4.3) 
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where 



Jr(si) OS JT(so) OS 



(4.4) 



a 



r(so,si) 



s dt 



dA 



so \Jr(s) 



1 



dA I ds. 



Remark 4.1. Using (2.9) - (2.6) we note that 



= tr (k'cos^ + r/tsin^) + n— . (4-5) 

We conclude that in the case of uncurved tube k, = t = 0, the right hand side 
of (4.3) vanishes, and the we get L{so,Si) = for all Sq,Si G (0,1). The 
equation L{sq, Si) = can be regarded as an integrated form of Webster's 
model. 



To obtain Webster's equation (3.4), we first study the hmit 



L{s,s') 



s'^s s' — S 



hm 



(4.6) 



here L{s,s') is given by (4.4). Recall that we aim at deriving (3.4) for the 
average 



1 



Ais) 



T{s) 



The first two terms in (4.4) give us the elliptic part of Webster's equation. 
We rewrite 

(j)dA] =^ {A{s)^) 



d_ 
ds 



as 



A{s) 



r(s) 
2R'(s] 



ds R{s 
= -A'(sU + 



ds 



(pdA + 



ris) 
d 



d_ 
ds 



(pdA 



ris) 



OS \Jris) J 



The last term is the limit as s' — ?■ s of 
1 



s — s 



(pdA - / (pdA 
r(s') Jr{s) 



2- 0(s',r,^) -0(s,r, I 



^0 



s — s 



+ 



R{s') - R{s) 



s' — s 



2-K 



R{s')-R{s) 



-r dr dO 

"Ris') 
Ris) 



fs, r, 9)r dr dO, 
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which gives 



9_ 
ds 



(pdA 



ds 



'r(s) / Jr{s) 
From these expressions, we now obtain 



2n 



dA + R{s)R'{s) / (j){s,R{s),e)de. 



j^^ ^ ^^dA = A'is) (^0 - i- <p{s, R{s), e)de^ + A{ 



ds 



Thus 



hm 



1 



d(j) 



s'^s s' -s Vir(s') ds Jr(s) 9s 



ds 7r(s) ds 



dA 



(4.7) 



2n 



{s,R{s),9) d9 



and the reader can recognise the elhptic term of Webster's equation on the 
last hne. 



For the boundary dissipation term in (|4.4|) we get 

d 



a 



r{s,s') ^ 9t 



-dA = —a 



dt 



£ (^J%(~s,R(~s),e)^^Wi~s)dS 



by using (2.8). Thus 
hm ^ 



1 dcp 



dA 



s'^s s' -s Jr{s,s') 2 dt 



—a hm 



s — s 



dt 



«W^(5)^^ 0(s,i?(^),^)^j. 



To proceed we note that by (2.6), the expression inside the parenthesis takes 
the form 

(f){s,R{s),e)— = / (j){s,R{s),e){l-ri{s)cos 
^ Jo 

2tt 



--2'K(t)(s) - 27r 



1 

2^ 



[s,R{s),e)de 



2lT 



r]{s) / (f){s,R{s), 6) cos Ode. 
'o 
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Thus we conclude that 



-a 



hm — 

s'^s s' 



+ 27raW 



1 90 

r{s,s') ^ 9t 



dA = -27iaW 



dt 



d_ 
dt 



+ aWrj 



d_ 
dt 



AS) 

2-K 



1 

[s,R{s) 



2-K 



(j){s,R{s),e)de 

cos 6d6 



We take the by now famihar hmit also in the final term in (4.4). We get 

1 d\ 



lim 

s'^s s' — S 
- 1^ 



(j)dA 



2-^2 



dA ds 



1 



r(.) 



dt' 



-dA 



(4.9) 



T{s) - 



In order to obtain (3.4), we must express most of the contribution of the term 
containing Jp^.^^ in (4.9) using the averaged solution 0. Unfortunately, 
the curvature factor S is not constant, and we cannot just bring it out from 
under the integral sign. To deal with this problem, we use the sound speed 
correction factor S(s). It is clear that S(s)~^ is the best estimate for function 
^ over r(s) in the sense of least squares. We have 



S(s 



|2 ■ 



A{s) 



V{s) 



dA 



(4.10) 



7ri?(s) 



2-n pR{s) 



^0 



[1 — rft(s) cos Oy rdrdO = 1 + 7?7^(s) 



where the curvature ratio ri[s) is given by (2.5). The error function in (3.14) 
satisfies the identity 



1 



With the aid of this we get the splitting 

(t)dA A{s) - 



V{s) 



+ 



E(l)dA 



V{s) 



which implies 



1^ 



T(s) - 



(j)dA\ _ A{s) d^(f) 
~ cisV dt^ 



r(s) 



Ed^ 



dA. 



(4.11) 
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Here 



c(s) := cS(s) 



(4.12) 



is the variable speed of sound. 

We consider next the terms on the right hand side of (4.3). For the terms 
containing the torsion r, we get 



hm 



s'->s S — S 



T(S 



r(s') 99 



dA - t(s) 



r{s) 



de 



dA 



r [s 



r(s) 

27r 



T{s) 



86 



dA = t{s] 



d_ 
ds 

d_ 
ds 



r{s) 



ris) 



d(j) 



dA 



(4.13) 



09 



dA 



because J^" %d6 = 0. Then we proceed as in the first two terms of (|4.4|) 
with II in place of 0, and obtain 



d_ 
ds 



ris) 



d(j) 



dA 



V{s) 



dsde 



dA + R{s)R'{s) 



2tt 



09 



{s,R{s),9) d9 = 0, 



^irnpp ^ - and f^'' ^ d9 - 

^1^^^ dsde - 89 ds Jo laJ - ^■ 

We conclude that only the first term remains on the right hand side of 



(4.3). Taking the limit, we get 
1 



lim 

s'->-s S' — S 



n{s,s') 



V(l)dV 



T{s) - 



VcpdA. (4.14) 



We now have everything we need for the limit (4.6). We put (4.7), (4.8), 
(4.9), and ( 4.11[ ) together, and obtain 



1 , L(s,s') 1 d 
lim 



A{s) s'~^s {s' - s) A{s) ds \Jris) 9s 

1 



dA 



lim 



a 



1 90 



dA 



1 



A{s) s'^s s' - s 7r(.,s') 2 dt A{s)c^ dt^ Vir(.) J 



(4.15) 

(f)dA\ 



1 d 
A{^d^ 



A{s) 



ds 



27raW{s) d<p 



1 



- F{s,t) -G{s,t) 



+ 



+ 



aW{s)i]{s) d 



A{s) dt V7o 



2n 



A{s) dt c2S(s)2 dt^ 



[s,R{s),9) cos 9d9 



A{s) 



EA^dA 



ris) 
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where F and G are given by (3.11) - (3.12). Equation (4.15), together with 



the wave equation (4.3) and equations (4.13) and (4.14), now gives Webster's 



equation (3.16). The final error term H in (3.16) is given by (3.13). If 



we assume that the right hand side of (3.16) is neghgible, we obtain the 



Webster's horn equation (3.4) with curvature and dissipation. 



Webster's control/observation boundary conditions at s = 



We next derive the boundary conditions (3.8) for Webster's equation at s = 



that correspond to the boundary conditions of the wave equation on r(0) in 



(3.6). 



We get for smooth by using (4.1) 



r(o) 
c 



cj^±=^]dA = ~c [ ^dA±^f <pdA (4.16) 



r(o) 



(0) 



The last term on the right hand side of (4.16) vanishes if k(0) = which 
will be assumed. Thus, our interest lies in the term 

dA = I I —{s,r,9)rdrd9 for s = 0. 



r(o) 



We get 



ds 



271 pRiO) 



'0 Jo 
ds 



JO 



— (s, r, 6)rdrd6 
as 



\J0 



R{0) 



{s,r,6)rdr + / (f){s,r,6)rdr \ dO 
Ir{s) 



d{A4>) d ( /•^(') 



ds ds 



is, r, 9)rdrd9 



'0 J R{0) 

For the latter term on the right hand side we get 

ds 



(s, r, 9)rdr ) d9 



2tt 



<P{s',r,9)-<p{s,r,9) 1 
hm / : rdr + iim 



R(0) 



s' — s 



/ / -^rdrd9 + R{s)R'{s 

'o Jr{o) JO 



2lT 



s'-5-s S' — S 



{s,R{s),9) d9 



Ris') 



R{s) 



(s, r, 9)rdr ) d9 
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and hence by setting s = above we obtain 



r(o) 



^dA = (0) - Rio)R'{o) r '/»(o, ^(0), e) de. 

OS as In 



We conclude from this and (4.16) that if k{0) = 0, we have for all t > 

d(j) d(j)^ 



1 



A(0) 



r(o) 



'^(9z/ ^ dt 



dA 



c^{t,0)±^{t,0)-K{t) (4.17) 



ds 



dt 



where 



A{0) 



27C 



27r 



(j)it,0,RiO),e)d9 . 



Now, if satisfies (3.6) then satisfies 



(4.18) 



-c(0)^(0,t) + ^(0,t)-ir(t) = 2 



ds 



dt 



pA{0) 



U{t) 



- c(0)^(0,t) - ^(0,t) -K(t) = 2^^y(t) 



(4.19) 



where u(t) := -^j^^ f-p^Q^u{r,t)dA. The assumption k(0) = is convenient 
here, too, since it implies c(0) = c where c(s) := cS(s) is the sound speed 
corrected by S(s) in ( 3.3[ ). 



Thus, for smooth we have (3.6) =^ (3.8) if fi;(0) = A'{0) = which we 



make a standing assumption in p^; see also [221 P- 1992]. 



Proof of Theorem 3.2 



In this section, we give the proofs based on the formal computations in the 
previous section. Some auxiliary results are needed first. 



Proposition 5.1. Assume that Q is a tubular domain described by (2.3) 
and : — 7- M satisfies the regularity assumptions cf) G C^(]R+; if ^(fi)) and 
A(p e C(]R+;L2(f2)) m (|L6|). ByV ddVl denote the walls of the tube. Then 
the boundary trace satisfies E C (R.'^ ; (T)) . 



Proof. Let ip G H^{Q) such that Aip G L^(fi). Since the radius function 
R{s) is smooth, we may assume that the tube Q is of constant diameter near 
its ends r(0) and r(l) — if not, use a diffeomorphims to obtain R'{s) = 
for s G {0,r]) U {1 — rj, 1) for rj > 0. We extend Q from both ends to a 
longer tubular open set Q that has a smooth boundary. This extension can 
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be carried out in many ways but the function ^jj must be extended to all of 



Q so that (1.6) are satisfied by the extended function ip. Let us consider the 
end s = 0. 

For s G {—ri/2, 1) define the extension by reflection 

^(.,r,^) = <('^^^'^'^^ for. G (0,1) 

^ ^ ^ 3^{-s,r,9) -2ij{-2s,r,9) for s G (-1/2, 0) ^ ^ 

and similarly at the other end. This extension gives us a function ip defined on 
a tube of length l + rj, and by smoothing the tube around the ends produces 
Q having a smooth boundary with f2 C and F C dQ. It is easy to see that 
ip satisfies ( 1.6[ ) since ip does. 



Thus ip G H\n) with Aip G ^^(fi) and g G L'^{dn) where the boundary 
dCl is smooth. By elliptic regularity theory (see Remark 7.2 on p. 188]) 
(with s = 3/2, m = l, uiq = 1), we conclude that ip G i7^/^(n). Using 
[H Theorem 7.4] we conclude that G i:'i^^(l^) := {/ G H^/^{n) : A/ G 
S~^/2(fi)} where has been defined on [HI p. 172] since L'^{^1) C 

S-V2(f]). It follows from [H Theorem 7.4 on p. 188] that i>\g^ e H^{dn) 
and thus V'lp ^ H^i^) by restriction. Now, the trace mapping h-)- ip^^ is 
bounded from E(A;L2(r])) := {/ G if^(fi) : A/ G ^^(fi)} into L2(fi), and 
its range is in if^(F); it is thus bounded from E{A] L'^{Q)) into H^{r). Since 
(j) G C(M+; E(A; L^jfi))), we conclude that 0L G C(M+; //^(F)). □ 



Proposition 5.2. Let F 6e as m Proposition 5.1. Define for g G H^(r) the 
linear mapping B by 



{Bg){s):=—\ gis,Ris),e)de for s G (0, 1) 







Then 

B & C{H^{T)-H^{Q,l)) for all A; G M. 
Proof. For any g G ^""(F) and w G /.^(0, 1), we have Bg G C°°(0, 1) and 



V^O 



(%,w^)l2(o,i) = ^ / (/ ^(s,i?(s),^)d^ ) ?i;(s)ds 



1 /.27r 



^0 



g{s,R{s), 



2ttW{.s 



where dA = W{s)d9ds and w{s,9) = 2wW{s) ■ Because C°°(F) is dense in 
H\T), we conclude that B* G £(if-'=(0, 1); i/~'=(F)) for all A; > where 
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{B*w){s,R{s),e) = 2^!^ for all w G ^^(O,!)) and {s,R{s),9) e T. In 
particular, B G C{Hl^{T); H^{0, 1)) for all A; > 0. 
Define the spaces 

W={re C{T) : r{s,R{s),e) = sro{e) + (1 - s)n{e) for r^n G C{T)} 

and FT'' := WnH''(T) equipped with the norm of H''(T). Then is a closed 
subspace of if''(r) for k >2, and it is easy to see that W^+Hq(T) = H''(T). 
The operator B maps W'' boundedly onto linear functions on (0, 1) that are 
equipped with the norm of H''{0, 1). We conclude that B G C{H''(T); H''{0, 1)) 
for k > 2 and k = 0. The claim follows from this by interpolation; see, e.g., 
p3l Theorem 5.1 on p. 27 and Theorem 7.7 on p. 36]. □ 



Proposition 5.3. Let Q be as in Proposition 5.1 For f G L'^{^), define the 
linear mapping A by 

{Af){s):=^[ <pdA for SG(0,1). 

Then 

A G c{L^{Si)] l2(o, 1)) n c{H\ny, h\o, i)). 

Proof. For G L'^{Q) we have by Holder's inequality 

< / ( ! <i^''dA \ ( f dA] ds 



r(s) J vr(.) 



^ ~'s) ,9 rdrdOds ^ 

^ ' ' " '^"L2(f7) 



where C\ := maXsg[o,i] because rfA = rdrdO on r(s) and dl^is given by 



(2.7). Using the operators A and B, equation (4.7) takes the form 

first for G C°°{U), and then by density for all G if^(fi) and G H^^^{T) 
and P Theorem 1.5.1.3]. It follows that ||^ M0)||l2(o,i) < C^illt • V0||l2(o,i) + 

sup,,[o,i] ^ (C2||0|rlU^(r) + C^i||0||l^(o,i)) < CUUhu) for all G i^i(fi) 
for some C < oo. □ 
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We are now ready to give a rigorous proof for Theorem 3^ 



Proof of Theorem \3.S\ Claim (|ij): By (1.6) we have 
1^ G C^(]R+; -L^(0, 1)). Thus an apphcation of Propositions 



5.3 



L^jn)) and 



and 



5.2 



im- 



phes that e C2(M+; ^^(0, 1)) and f G Ci(M+; ^^(0, 1)), as desired. 

Claim The functions A(s), A(s)-\ 77(5), Vr(s), S-^V(H-^), a nd ^ 
are all smooth by assumptions. Hence, by inspection of formulae (3.11) 
- (3.13) we need to show, in addition to claim (pi) of this theorem, that 



as Wr) e C{1 



L2(0, 1)) in order to prove that F G C(M+; L2(0, 1)) and 



G,H e Ci(M+; L2(0, 1)). This follows directly from Propositions 5.1 and 5.2 



Claim (iii): By 1.6 and (3.15), all the computations in Section |4| hold 
pointwise almost everywhere, apart from (4.7); might not have two deriva- 
tives with respect to s. To deal with this, we interpret the second derivative 
in (4.7) in the sense of distributions]^ 

For |/;,| < s < 1 — we may write, in the notation of section |4| 



L(s, s + h) 



■ V(pdV, 



(5.3) 



where 



L(s, s + h) 



T{s+h) 



d(f) 
ds 



dA 



V{s) 



d(f) 
ds 



dA 



— a 



dt 



dA- 



1 



dt 



-dAdo. 



'r{s+h,s) 

Note that here we used the assumption that the torsion r vanishes. Take now 
an arbitrary function ( G C^{{0, 1) x (0, T)), and choose h sufficiently small, 
so that the spatial support of ( is contained in the interval (2|/i|,l — 2|/;,|). 
By a simple change of variables, we see that 



r{s+h) 



ds 



dA({s,t)ds 



r{s) 



d(j) 
ds 



dAr]{s — h)ds. 



(5.4) 



Next we multiply (5.3) by C, integrate over s, divide by h and use (5.4). This 
leads to 



^ L{s,s + h) 



({s, t)ds 



— a 



V{s) 

1 
h 



ds 



h,t)-as,t) 

h 



ds 



r{s+h,s) 
s+h 



r{a) 



1^ 

Edi 

1 d^i 

'^di^ 



dA({s, t)ds 



dAdads. 



^The (iin)patient (reader) may now skip the rest of the proof... 
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Then let h tend to zero. The hmit of the first term on the right is 



hm 

/i-s>0 



Vis) 



OS h 



r(s) 



OS OS 



In the other terms, we may pass to the limit h —> (or, equivalently, s' — )■ 
s) as in Section 111 since all the formulae after (4.7) hold pointwise almost 



everywhere. An integration in time then leads to (3.17). 



Claim (iv): As argued in Section |4| we have (3.6) =^ ( 3.8[ ) for smooth 
if KiO) = R'iO) = 0. Since §J(-,t) e L^m) and^|^(^^(-, t) G H'/'iTiO)) 



by (1.6), and the claim follows because the averaging over r(0) is bounded 
linear functional on L^(r(0)). 



□ 
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A The surface area element 



By (2.3), the wall V can be parametrized as 

r (s, 61) = 7(s) + R{s) (n cos 6^ + b sin 6) 

where s G [0, 1] and 9 G [0, 27r]. Hence the surface area element can be found 
by using the formula 

dS=\f^yi g| dsdO. 



We have |^ = -R(s)(— n(s) sin6' + b(s) cos6'), and by (2.1) we get 
Qr 

— = t(l - kR) + n(R' cos 6 - Rt sine) + h( R' sine + Rt cos 6). 

OS 

Since {t, n, b} is a right hand orthogonal system, we see that 
X n = b(l - kR) - t( R' sine + Rt cose) 

OS 

and 

X h = - n(l - kR) + t( R' cose - Rt sine). 

OS 

Hence the exterior normal derivative lies in the direction of 

- X = R\-tR' + n(l - kR) cose + b(l - kR) sine] (A.ll 

OS de 



since kR < 1. Now (2.8) follows by orthonormality, see (2.5) 
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B The volume element 



To verify (2.7) we need to compute tlie Jacobian of tlie coordinate trasfor- 
mation: 



d{s, r, 6) 
d{x,y,z) 



ds 
dx 
dr 
dx 



ds 
dy 
or 
dy 

m 

dy 



ds 
dz 
dr 
dz 

m 

dz 



[-rSt 

I COS 6 



dx 

H(t ■ i) 

[cos n + sin b] ■ i 



s(t-j) 

[cos ^ n + sin 6* b] ■ j 



+ S sin ( 



S cos^ 9 



sin 6 _j_ cos 6 y^l 

r r } 


i 


t ■ i 






n ■ i 






-rSt - -\ 


cos 6 
r 


t i 






bi 






-rSt - + 

r 


cos 8 1 
r 


t • i t ■ j t • 


k 




n ■ i n ■ j n ■ 


k 




b ■ i b ■ j b ■ 


k 





-rHt - 



-n 



tj 
nj 

b] ■ i [-rSt - + 

tj 
bj 

□1 ■ i T-rSt - + 



S sin^ e 



t ■ i t • j t • k 
b ■ i b ■ j b ■ k 
n ■ i n ■ j n ■ k 



S(t • k) 
cos ^ n + sin b] ■ k 
^b] -j [-rSt - + e^^b] 

tk 
n ■ k 

^b]-j [-rEt 

tk 
bk 

b]-j [-rSt-^n- 

t • i t ■ j t ■ k 
n ■ i n ■ j n ■ k 
b ■ i b ■ j b ■ k 



cos t 
r 



sin 9 ^ I cos 9 

r r 



cos 9 



b] 



We write tlie basis change as xi + yj 
that 



zk = at + (3n + 7b, and thus we see 



a 




X 




't ■ i 


t 


j 


t 


k 


/3 


= A 


y 


where A := 


n ■ i 


n 


■ j 


n 


k 


7 




z 




bi 


b 


• j 


b 


k 



Because both the bases {i, j, k} and {t, n, b} are orthogonal, the basis change 
matrix A is unitary. Hence its determinant is of absolute value 1, and 

d{s,r,9) ' 
d{x,y,z) 



-. From this we conclude that (2.7) holds 



C The gradient and the normal derivative on 

r 

We have 

„ . d . d ,9 

^ = ^7r +J7^ + ^^7^• 
o'x ay oz 



Suppose that (j) = (j){s, r, 9) where s, r, and 6 are functions of x, ?/, and z. 

d 9i_d_ _|_ dji _d_ I d9 d 

dx dx ds dx dr dx d9 ' 



We have by the chain rule -4- = f^^ + l^^ + l^i^, and hence we need to 
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compute 1^, 1^, and 

dX-—,OX ' ox 

by using (2.1 ) 



By differentiating (2.2) with respect to x, we get 

„x ds 



[L — rn cos t/ 1 — t(s) 
ax 

-rr sm c^t; — h cos o— r sm n[s) 

ox ox J 



dx 



+ rrcosu- — h sm&— — h rcosu— b(s). 
\ Ox Ox Ox J 

Because t(s), n(s), and b(s) are orthonormal, we get 



1 — TK, COS 6 










- ds - 

dx 




'tis)- 


—rr sm6 


cos 6* 


— sin^ 




Or 




n(s) 


rr cos 6 


sin 6 


cos 6* 




r% 

. dx. 




b(.)_ 



By the topmost row, |^ — -, , 

-t^ ' ax 1— rrecost 



t(s) 



i, and using this gives 



dr 

dx. 



COS 6 sin 6 
- sin 9 cos 9 



n(s) + 



h{s) - ^ 



-Ms) 



1 — TK COS C 

rr cose 



rK cos S 



or 



t(^) 



■ 1 



l—rn cos ( 

[cos 9 n(s) + sin 9 b(s)] ■ i 

^ . t(g) - ^n(s) + 

L 1— rKcos9 ^ ' r \ / ' 



'-Us 



We now conclude that ^ = F ■ i where 

dx 



F{s,r, 




+ n s 



A similar argument shows that also 



dy 



F-j and 



F ■ k hold, and hence 



V0 = -F. The formulas (2.9) and (2.10) can be read from this. 

It remains to compute the normal derivative on the tube wall T. By 



dAll ), (Q, and (Q we get 

W 
R 



W 



tR' + KR){n cos + b sin 9)) ■ {tDi + nD2 + bDg) 

d 



-R! 



d__ d_ 
ds~^d9 



+ (1 - kR) 



dr 



as desired. 
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D The Laplacian 

Let us proceed to compute A = in tube coordinates. We have 

= (tDi + nD2 + bDg) ■ {tDi + nD2 + hD^) 
= t ■ (DitDi) + t ■ (Z^inDa) + t ■ (/^ibDg) 
+ n ■ (Dat/^i) + n ■ (DsnZ^s) + n ■ (DsbDg) 
+ b ■ {D;tDi) + b ■ (D3n/^2) + b ■ (Z^sbZ^a). 



(D.l) 



The first term in (D.l) takes the form 



t ■ {DitDi) = t ■ [iDit)Di + tDl] = Dl 



since (1 — rK,cos9)Dit 



dt 



de 



n — = n which is orthogonal to t. 



The second term in (D.l) gives 



t ■ (DinDa) = t ■ [{D^rv)D2 + rvD^D^] = t ■ ((1^111)^2) 



TK cos 9 



because (1 — rKcos^^)-Din 



dn dn 

ds de 



-nt + rb and t ± b. 



The fifth term in (D.l) is as follows: 



n ■ (DznDa) = n ■ [(Z^2n)Z^2 + nD^] = 
since An = cos - = by virtue of the fact that n = n(s) 



The ninth term in (D.l) goes like this: 

b ■ (D^hD^) = h ■ [(A3b)I^3 + 



^3 



since D^h = sin 6*^ + = by virtue of the fact that b 



Us). 



All the remaining terms in (D.l) vanish by the facts that = D2n = 
= D^t = -Dsn = D^h = and orthogonality of the vectors t, n, and b. 
The laplacian now takes the form 



Dl + Dj + Dj- kED 



(D.2) 



where the operators Di, D2, and are given by (2.10), the curvature is 
denoted by k, and the curvature factor S is given by (2.6). We now have by 
(|2l0l) 



(9s 2 



2r- 



92 



92 



dsde ' ' 06^ 

d ^ d 
rE(KT sin 6 + k' cos 6)— — h (rrEiKT sin 9 + k' cos 6) — t') — ; 
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D'^-cos^e^'^ sin 2^ d"^ sin^ ^ 9^ sm^ 6 d sin 2^ 9 
2 Q^2 J, QrdO dO"^ r dr 89^ 



and 



= sin^ e 



sin 2^ 92 cos^^ cos^ 9 sin 20 9 



+ 



+ 



r2 



+ 



r 9r 



r2 90 



which imphes 



^2 + Dl 



92 1 92 19 

Qy,2 j,2 QQ2 J, Qj, 



Thus by (D.2) 



92 



92 



A = ^ h 

9s2 9r2 



9 



rS (/tr sin + cos 0)^j — h I /tScos0 



9s 



+ ( rT'E?{KT sin 9 + k' cos ( 



52 „2 92 

2rS^ 

902 9s90 

9^ 

9r 

KSsin0\ 9 
90' 



(D.3) 



We remark that if 7 is a plane curve, then r = 0. If, in addition, k = 
(contrary to our standing assumption) we get S = 1, and then (D.3) becomes 
the laplacian in cyhnder coordinates 



92 92 19 1 92 
A = — + — + -— + 



92 19/9 
+ -— I r- 



Qg2 Qy,2 Qy, y,2 QQ2 Qg2 r dr \ dr 



+ 



1^ 

^902 



We conclude that formula (D.3) is valid even if k = 0, even though not all of 
the Frenet formulas make sense. 
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